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Abstract 



Power-suppressed corrections to parton evolution may affect the theoretical 
accuracy of current determinations of parton distributions. We study the role 
of multigluon-exchange terms in the extraction of the gluon distribution for 
the Large Hadron Collider (LHC). Working in the high-energy approximation, 
we analyze multi-gluon contributions in powers of 1/Q 2 . We find a moderate, 
negative correction to the structure function's derivative di^/dhiC} 2 , char- 
acterized by a slow fall-off in the region of low to medium Q 2 relevant for 
determinations of the gluon at small momentum fractions. 
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The estimate of the theoretical accuracy on the determination of parton distributions is 
relevant for phenomenology at the Large Hadron Collider (LHC) [1-3]. A potential source 
of theory uncertainty is given by power-suppressed contributions to parton evolution. In 
particular, the extraction of the gluon distribution for small momentum fractions depends 
on data in the region of high energies and moderate to low Q 2 , in which power-like corrections 
from multiple parton scatterings are potentially significant. 

The purpose of this note is to derive an estimate of these corrections, based on high- 
energy amplitudes for multigluon exchange [4]. We concentrate on the power correction A 
to the Q 2 derivative of the structure function F 2 



where K is the perturbative kernel, G is the gluon distribution, and A is 0(1/ Q 2 ). This 
correction contributes to the theory uncertainty on G, since below x ^ 10~ 2 this is mostly 
determined from data for dF 2 / dhiQ 2 . 

The theoretical framework to treat multiple scatterings is based on the s-channel picture 
of DIS [4] , and its basic degrees of freedom are described by matrix elements of eikonal lines. 
The corresponding predictions, incorporating nonperturbative dynamics, are valid down to 
small Q 2 . To enforce consistency with the framework of standard parton analyses [2,3,5] we 
will expand the answer in powers of 1/Q 2 - We study the behavior of the power expansion 
at low Q 2 and x, and identify the power corrections by subtraction of the leading-power 
contribution. 

We find moderate, but non-negligible negative corrections, increasing in size as x de- 
creases. We find that with a physically natural choice of parameters in the eikonal matrix 
elements we can achieve a sensible description of data and still have power corrections to 
dF 2 /d\nQ 2 that do not exceed the leading power already at Q 2 ^ 0.5 GeV 2 . However we 
find that for small x the corrections have a slow fall-off with Q 2 in the region of intermediate 
Q 2 , Q 2 ~ 1 -h 10 GeV 2 , behaving effectively like l/Q x in this region, with A close to 1 for 
x ^ 10~ 3 . This behavior results from summing the power expansion. As a consequence the 
power corrections stay larger than 10% up to Q 2 of a few GeV 2 for x below 10~ 3 . 

The contents of the paper is as follows. We first give the s-channel results that provide 
the basic elements for the evaluation of the power expansion. Then we focus on the Q 2 
derivative of the structure function. We examine the next-to-leading-power contribution 
and the sum of the power series, and present numerical results. 

We begin by recasting the result [4] for multi-gluon contributions to DIS structure func- 
tions in the Mellin-transform representation. This is convenient to analyze the 1/Q 2 expan- 
sion. In the high-energy approximation we describe gluon exchange in terms of eikonal-line 
operators 



where A is the color potential and V is the path-ordering. Following [6] we notate the matrix 
elements of eikonal operators at transverse positions b and b + z as 



dF 2 



K®G [1 + A] + quark term , 



(1) 
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S(z,b)=|[dP'] (P'\^-Tr{l-V\b + z)V(b)}\P) . (3) 

Here [dP'} = dP' + d 2 P' /2P' + (2n) 3 , z is the transverse separation between the eikonal lines, 
and b is the impact parameter. The result [4] for the transverse structure function F T can 
be written as 

pc+ico An, f _ 

xF T (x,Q 2 ) = P: / d 2 h E(u,b) , (4) 



J c— too Z7T2 ^ 

where < c < 1, S is the Mellin transform of H 

2Kb)=/0(z 2 r 1 S(z,b) , (5) 
and $ is a calculable coefficient. To lowest order 

*M -Vc 2 ^ i (c?2r ^ 2 r(3- M )r(2- M )r(i- M )r(2 + M ) 

W V a 16 7 r% Wj 4« r(5/2- M )r(3/2 + M ) ' l j 

where N c = 3, and e a is the electric charge of quark of type a. 

From this approach the parton-model framework is recovered through an expansion in 
powers of gA, valid for small z. In particular, the coefficient of the quadratic term in the 
expansion of E can be related to the gluon distribution G [4], 



d 2 b E ~ (xG) z 2 (1 + 0(\z\)) , (7) 



AN. 
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where T R = 1/2. 

The quark distribution can be dealt with by the same method. The main difference 
is that while in the structure function case the ultraviolet region is naturally regulated by 
the physical scale Q 2 , in the case of the quark distribution we need to treat the ultraviolet 
divergences. The result in dimensional regularization is [7] 

xq(x, n) = Gu 2 )~ 2e J d 2 - 2€ z d 2 - 2e b w(z) S(z, b) - UV , (8) 

where 

, , N c 1 /^ 2 z 2 \ 2e r(2-e) 2 

and we have indicated by UV the ultraviolet subtraction. This is required since in Eq. (3) 
(1 - V*V) -> for z -> and S oc z 2 . Using Eqs. (7), (9) the x < 1 form of the evolution 
equation for the quark is reobtained from Eq. (8). Power-suppressed contributions arise 
from the difference between the terms left over from the ultraviolet regularization in the 
quark-distribution and structure-function case [7]. In general, these depend on the scheme 
used for the ultraviolet subtraction. This dependence does not enter in the correction to the 
structure function's derivative which we consider next. 
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The matrix element S is nonperturbative and is to be determined from experiment. 
For the numerical calculations that follow we model its functional form according to the 
model [8,9] 

*-">-^(^r ' (10) 

where /x s is the saturation scale, with b dependence as in [8]. The operator relation (7) 
implies for model (10) that 



/ 



dh $ = a *M xG{x,iif) . (11) 



c 



In Eq. (11) we have indicated explicitly the dependence of the running coupling and gluon 
distribution on the renormalization/factorization scales /i r , fj,f. In the present context the 
choice of these scales amounts to specifying the model for S. 

Consider now the derivative of Ft with respect to InQ 2 , F' T = dFx/ dhiQ 2 . Taking the 
derivative cancels the factor Xju in Eq. (6). We determine the expansion of F' T in powers 
of 1/Q 2 by closing the integration contour in the complex -u-plane to the left and evaluating 
the residues at the poles of the integrand. First we verify that the result from the leading 
pole (LP) u — coincides with that from Eq. (8) for the quark distribution. We get 

N r 

^t,lp = E J d " h Res «=o 3 • (12) 

By inserting Eqs. (10), (11), Eq. (12) yields the perturbative leading-power coefficient. We 
identify the power-suppressed correction by subtracting off this contribution. Next we con- 
sider the contribution from the next-to-leading pole (NLP) u = —1, F^ NLP . This is propor- 
tional to the u — —l residue of S. In Fig. 1 we compute the ratio 

£(nlp) = Fy jNLP / (F ThP + Fy jNLP ) (13) 

versus Q 2 at x = 10 -2 and x = 10~ 4 for two different choices of /i r , fj,f. The natural scale for 
fi r and \Xj should be set by the inverse of the mean transverse distance z. For the illustration 
in Fig. 1 we take this scale to be on the order of Q, and plot results for fif = Q, fi r = Q and 
jif = 2Q, fi r = Q/2. We use the CTEQ parton distributions [3]. 

The ratio #( NLP ) goes like 1/Q 2 - We see from Fig. 1 that the size of the NLP contribution 
is rather sensitive to the scales fif and fi r . This is not surprising. The definition of the model 
for S, embodied in the choice of fif, fi r , corresponds to defining the (otherwise arbitrary) 
separation of perturbative and nonperturbative effects. In particular, varying these scales 
amounts to effectively simulating contributions of higher perturbative order. The variation 
of the nonperturbative power correction in Fig. 1 says that this is unambiguously defined 
only once we specify what we include in the perturbative part of the calculation. See [10] for 
an extensive study of the issue of scale-setting and running coupling effects in high-energy 
evolution. In what follows we simply set the scales from comparison with experimental data. 

Beyond the next-to-leading power, the poles in the w-plane have multiplicity higher than 
1, leading to InQ 2 enhancements of the power corrections. The order-n term is of the form 
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FIG. 1. Relative size of the next-to-leading- and leading-pole contributions to the Q 2 derivative 
of the transverse structure function. 

c{nMQ2) Wf ' (14) 

where £ n give the dimensionful nonperturbative scale in terms of the b-integral of the mo- 
ments (5) of S, and C are coefficients determined from Eq. (6). Through the moments of 
S the correction (14) receives contribution from the exchange of any number of gluons via 
eikonal operators. 

We now proceed to evaluate numerically the contribution of all powers in 1/Q 2 . We 
go to z-space via the inverse Mellin transform and let the scales fif and fi r vary with the 
distance z. Details for the numerical integration in coordinate space are given in [7]. Adding 
to Eq. (4) the contribution of longitudinal F L , we compute the structure function F 2 and 
tune the factorization/renormalization scales, /if = Ci/\z\ and /i r = c 2 /|z|, by comparing 
the answer with the experimental data [11]. The result of doing this is reported in Fig. 2, 
where C\ = 4, c 2 = 0.32. 

Using these values for the model parameters, we turn to the derivative dF 2 /d\nQ 2 . We 
calculate the power correction by subtraction of the leading power as described around 
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FIG. 2. Predictions for the structure function F 2 compared with the data [11]. 

Eq. (12). In Fig. 3 we plot the result for the power correction normalized to the full answer 
and multiplied by (—1). We see from Fig. 3 that above Q 2 = 1 GeV 2 the corrections are 
below 20% for x ^ 10~ 4 . We take this as an indication that the power expansion is not 
breaking down, at least to such values of x. In fact, we find that corrections do not exceed 
the leading power already at Q 2 ~ 0.5 GeV 2 for x above 10~ 5 . 

Fig. 3 also shows, however, that as x decreases the power corrections remain non- 
negligible up to higher and higher Q 2 . For x < 10~ 3 it takes Q 2 of a few GeV 2 before 
the correction is less than 10%. This is associated with the curves having a rather slow 
decrease in the range of medium Q 2 in the figure, much slower than the asymptotic 1/Q 2 - 
For instance, for x ~ 10~ 3 the behavior in the region Q 2 ~ 1 4- 10 GeV 2 is closer to an 
effective power l/Q x with A ~ 1.2. This slow fall-off results from summing the terms (14). 
We leave to a future study the question of whether this behavior can be interpreted in terms 
of an effective, x-dependent semihard scale on the order of the GeV. 

In summary, the results above indicate that the power expansion should still work at 
x and Q 2 as low as in the region of the data presently used for determinations of the 
gluon distribution G for the LHC, but subleading corrections to dF 2 /d\nQ 2 from multigluon 
exchange are non-negligible in this region and contribute to the theory uncertainty on G. We 
have also computed the analogous corrections for the derivative of the transverse structure 
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FIG. 3. Power corrections to d^/dlnQ 2 versus Q 2 at different values of x. 



function [7], and we find that these are generally smaller than in the case of F 2 . This may 
be regarded as yet another motivation for the importance of a separate measurement of the 
longitudinal component F L [12]. 

A word of caution is needed in interpreting these results. Multigluon amplitudes are 
treated in the high-energy approximation. Also, the modeling of the nonperturbative matrix 
elements and the summation of the power series expansion call for a firmer understanding. 
Besides, power corrections from sources other than that considered here may be relevant as 
well. In particular, corrections from self-energy graph insertions are still largely unexplored 
for flavor-singlet observables [13]. 

Nevertheless, the method presented above allows one to obtain an estimate of multi-gluon 
corrections which can be made consistently with perturbative evolution order by order. It is 
based on subtraction of the leading pole in Mellin space. This serves to specify the definition 
of the power correction. In this work we have been concerned with the contribution to the 
Q 2 derivative of the structure function, relevant for the extraction of G. But the approach 
can in principle be extended to evaluate corrections to F 2 itself, and to processes directly 
coupled to gluons. The latter will be especially interesting for studying multiple-scattering 
effects in the production of jet final states. 
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